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NOMENCLATURE 


A : Heat transfer surface area 
App : Bump base area 

Apc : Bump cross sectional area 
Ay, : Bump surface area 

Ar: Passage flow area 

A. : Plate surface area 

b : Bump base diameter 

C, : Cold fluid capacity rate 

Cy : Hot fluid capacity rate 

C, : Capacity rate ratio 

Cp : Specific heat at constant pressure 
D : Equivalent diameter 

Г: Басов ан но 

G : Mass velocity 

h : Heat transfer coefficient 

hy : Bump height 

1: Colburn factor 

k : Thermal conductivity 

Кы: Contraction loss coefficient 
K, : Expansion loss coefficient 
L : Plate length 

m : Fluid mass flow rate 

N : Bump numbers 

Ntu 


Nu : Nusselt number 


: Number of transfer units 


P : Wetted perimeter 

Pr : Prandtl number 

Q : Heat tranfer from hot fluid to cold fluid 
q Мелих 

Ке: Reynolds number 

Ey : Wall resistance 

5: Arc length of the bump 


St : Stanton number 


[kg/m?s] 
[W/m*K] 


[mj 


[W/mK] 


[W] 
[УУ / m?] 


[W/m*K] 


Í lol ñuid temperature 

t : Cold fluid temperature 

Ty, : Hot fluid bulk temperature 

tc : Cold fluid bulk temperature 

T : Hot fluid average temperature 
t: Cold fluid average temperature 
I, : Plate surface temperature 

tp : Passage height 

U : Overall heat transfer coefficient 
V : Fluid velocity 

v : Specific volume 

w : Plate width 

AT : Temperature drop 

p : Fluid density 

д: Material thickness 

v : Kinematic viscosity 

tt: Dynamic viscosity 

I, : Dynamic viscosity at wall temperature 
©. Exchanger effectiveness 


ç: Ratio of free flow area to frontal area 


Subscripts : 


c : cold fluid 
fl. hot fluid 
1 : inlet 


2: outlet 
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I. INTRODUCTION 


A. BACKGROUND 

In a heat exchanger, one of the important parameters controlling the net heat 
transfer from the hot fluid to the cold fluid is the surface area separating the two fluids 
across which the heat transfer takes place. When large quantities of heat are to be 
transferred, computations usually indicate a requirement for large heat transfer surface 
areas. Increasing the area necessarily means either increasing the total path length 
travelled by the fluids in the heat exchanger or simultaneously decreasing the diameter 
of tubes and increasing the number of tubes. 

Where high volume is not a limiting constraint, flat plates can be employed. 
However, use of flat plates to obtain the necessary surface may involve other heat 
transfer considerations in the device. 

In order to enhance the heat transfer characteristics of flat plate exchangers, 
some type of obstruction may be placed between the plates in order to cause the fluids 
to take a sinuous path between them. These obstructions are often in the form of a 
honeycomb structure but they inherently possess several problems. First, there 1s 
usually a straight path through the heat exchanger in the direction of flow thereby 
defeating the purpose of the obstructions. In addition, the sharp bends and corners in 
- the honeycomb structure can create undesirable dead spots within the heat exchanger. 
They may also represent weak points and areas of deleterious stress concentration in 
the structure which, upon heating and cooling, will often crack. 

Interest in heat exchanger surfaces with a high ratio of heat transfer area to core 
volume is increasing at an accelerated pace. The primary reasons for the use of these 
more compact surfaces is that smaller, lighter-weight, and lower cost heat exchangers 
are the result. These gains are brought about both by the direct geometric advantage 
of higher ‘area density’ and also a higher heat transfer coefficient for the smaller flow 


passages. 


ЕП THESIS OBJECTIVES 
e ms eS 15 tO determine the effect of protuberancts in a 


rectangular passage on the convective heat transfer coefficient. 


[n order to accomplish this objective, measureinents of fluid mass flow-rate, inlet 
and outlet temperatures and exchanger pressure drops were obtained for several 
different fluids flowing in a cross-flow compact type heat exchanger containing the 
protuberances. Heat transfer coefficients were then obtained. From these data, 
correlations for Nusselt number, Colburn factor and friction factor were then 


developed. 


H. THEORY OF OPERATION 


A. GENERAL 
The rate of heat transfer between a solid surface and a fluid can be expressed in 


much the same way as for simple heat conduction through solids: 
Q=hAAT (em 21) 


Here AT is the temperature drop between the bulk fluid and the wall, Æ is the heat 
transfer surface area and fA 1s the heat transfer coefficient. The last term is equivalent 
to the term k/L in the simple heat conduction equation. 

Equation (2.1), which 1s often referred to as "Newton's Law of Cooling” is casy to 
use except for the evaluation of the heat transfer coefficient A, which depends on the 
fluid flow conditions, the thermal properties of the fluid, and the passage size. 

[t is often convenient in heat exchanger work to combine the two surface heat 
tramsier coellieients (that ıs, for the hot and cold fluid streams) and the thermal 
resistance of the tube wall to give a single parameter: the overall heat transfer 


coefficient, U, which may be defined as 


U=Q/AAT (Gam. 222) 


B. INDIVIDUAL AND OVERALL HEAT TRANSFER COEFFICIENTS 


Equation (2.2) may be written in another form : 
E- UTAL) (Са 9) 


Pere q is the heat flux from stream-1 to stream-2, having bulk temperatures T, and T, 
ИСО Се Оче сас wansler coeficient. 

It is now time to consider the questions of what precisely is meant by the 
interface area, and what is the relation of U to the individual resistances to heat 


transfer represented by the solid wall separating the two fluids, the deposits of dirt that 


E 


may cling to the surfaces ofthe wall, and the low-velocity fluid films that adjoin those 
deposits or the wall itself. 

Fig. 2.1 illustrates the temperature profile discontinuities at locations a (inner 
surface of inner deposit), b (outer surface of inner deposit and inner surface of tube 
Wall), c (outer surface of tube wall and inner surface of outer deposit), d (outer surface 
of outer deposit). It can also be seen that the fluid temperatures lie both above and 
below T', and T,, for these are defined as flow averages, that 1s, averages for the duct 
cross sections in which each element of area 1s weighted by the specific heat C, times 
the local mass flow rate p V. 

Thus, 


Tior? * (J & C; V T 0s/ J p Cg V 05) оз (eqn 2.4) 


where the averages are taken over the whole of the cross section of the ducts in 


question. 





Figure 2.1 Temperature distribution in the neighborhood of a tube wall. 


Let A, ¿ and h > now represent the film heat transfer coefficient on the inside and 


outside of the tube, respectively; let h, and A. represent heat transfer coefficients for 
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the two dirt deposits; and A, represents the heat transfer coefficient for the tube wall 


SON un TR and М will all be related to the relevant material thickness ó and 


thermal conductivities k by 


hab = Kab / Sab (eqn 2.5) 
hpc = Ке / О с (ғап 226) 
(ean 2.) 


hed - Ked / беа 


Each heat transfer coefficient, to be meaningful, must be associated with an area 


A. When this is done, the fact that the heat flow through all the elements in question 


is the same leads to 


Beer Е. Аа, l5) (сап 2.8) 
= hab Aab (Ta ТЫ 

Abe Abe (Th — Tg 

= hed Аса (Те 7 Tg) 

= Rp Ag (Ta T T2) 


Here A), Abe Ag Can most appropriately be taken as the arithmetic means of 
the inner and outer surface areas of the materials in question, and A, , and A „ arc 
best taken as the surface areas pertaining to locations a and d respectively. The arca 4 
without subscript, associated with U, requires arbitrary definition; it can be set equal to 
any of the other areas. 

Elimination 7,, Т,, T., and Ty now leads to the following important connection 


between the overall coefficient U and the individual coeflicients 


ОИЕ Е ИИТ) E CAT eas) (A he Ad) (can, 2-9) 
ПРИ (A/hpAq) 


Evidently, the relationship is best understood as implying that the reciprocal of 
UA is the overall resistance to heat transfer, and that the overall resistance 1s equal to 
NR Tu — ] — | M 
the sum of the individual resistances (A, ,4,,) > (habhab) and men 
In the case of heat transfer across a plane wall, all the areas take ıdentical values. 


Then Eaton ZIFredugessis 


U=[(1/h,,) + (I/hi) t (1 бре + (17 а] ! (сап 2.10) 


C. HEAT TRANSFER COEFFICIENT DEPENDENCES | 

A dimensional analysis easity shows that the heat transfer coefficient may be 
found to be related to two dimensionless parameters : the Nusselt and Prandtl 
numbers. The Nusselt number, Nu, is the term used for the grouping AL/Kk, where L is a 
typical length dimension that pertains to the apparatus, and & is the thermal 
conductivity of the fluid. This quantity is the ratio of the heat transfer coefficient to the 
thermal conductance implied by the quantity k/L. It can be seen intuitively, from the 
simple conduction equation, that the heat flow rate to a fluid flowing through a 
passage should be proportional to the thermal conductivity divided by a representative 
distance in the direction of heat flow such as the passage equivalent diameter. 

Obviously, further prescription is needed of precisely what dimension is in 
question (such as duct diameter), and just where the thermal conductivity is to be 
measured (such as at the interface, at the bulk condition, or in some reference state). 

The Prandtl number, Praise the term Пс ер ето СЕ This quantity is the 
ratio of the molecular diffusivity of momentum (as indicated by the viscosity) to the 
molecular diffusivity of heat (as indicated by the ratio of the thermal conductivity to 


the specific heat). 


D. HEAT TRANSFER COEFFICIENT LIMITATIONS 

The concept of the heat transfer coefficient is so convenient, and so reminiscent 
of electrical circuit theory, that there is some danger that it’s validity may be presumed 
to be unlimited. A few remarks on it’s limitations are therefore in order. 

First, the heat flux (to or from the fluids) may not be proportional tome 
temperature difference alone. Other examples can be found, for example, Hes 


involving chemical reaction rather than phase change. 
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Second, the value of the coefficient is by no mean independent of the values of 
the property differences that multiply it. In the case of heat transfer by natural 
convection, it 1s "well known that the coefficient increases with the imposed 
temperature difference. 

Third, there are circumstances in which nonuniformities of temperature or 
concentration along the interface are so steep as to aflect the value of the transfer 
coefficient. In these circumstances, it may be best to work with the fluxes and 


property differences directly, avoiding all reference to coellicients. 


III. THE EMPRICIAL CORRELATION FOR THE CONVECTIVE HEAT 
TRANSFER COEBBIGIEDI 


The manner in which a convection heat transfer coefficient may be obtained 
experimentally is illustrated in Fig. 3.1. In a prescribed geometry, such as a flat plate in 
parallel flow, convection heat transfer occurs from the plate to the fluid, if the surface 


of the plate is at a higher temperature than that of the fluid. 


Е ПИ 





Figure 3.] Experiment for measuring the average convection 
heat transfer coefficient. 


The convection heat transfer coefficient A, which 1s an average associated with 


the entire plate, can be computed from Newton's law of cooling, 


О = р (eqn > 
where 
h : average convective coefficient [W/m^K] 
A. : plate surface area [012] 
T, : plate surface temperature [C] 
Tp : fluid bulk temperature [С] 


From a knowledge of the characteristic lengtlı Z and the fluid properties, the Nusselt, 
Reynolds, and Prandtl numbers can be computed from their definitions. The Nusselt 


number is a function of the Reynolds and Prandtl numbers. 


Nu = (hL/k) = К Ке, Рг) (eqn, 3:2) 
and 
Ве = VL/v | ешп?) 
Рг = C n/k (eqn 3.4) 
where 
L : characteristic length [na] 
k : thermal conductivity [W/mk] 
V : fluid velocity [m/s] 
it: fluid dynamic viscosity [Ns/m?] 
Co : specific heat [J/kgK] 
v : fluid kinematic viscosity [m/s] 


The forgoing procedure is reproducible for a variety of test conditions. One can vary 
the velocity V and the plate length L, as well as the nature of the fluid, using for 
example air, water and glycol, which have substantially different Prandtl numbers. We 
would then be left with many different values of the Nusselt numbers corresponding to 
a wide range of Reynolds and Prandtl numbers and the results could be plotted on 
log-log scale, as shown in Fig. 3.2. 

Each symbol represents a unique set of test conditions. As is often case, the 
results associated with a given fluid, and hence a fixed Prandtl number, fall close to a 
straight linc, which may be represented by an algebraic expression of the form 
Since the values of C, m and n are independent of the nature of the fluid, the family of 
straight lines corresponding to different Prandtl numbers can be collapsed to a single 


line by plotting the results in terms of the ratio, Nu/ Pr", as shown in Fig. 3.3. 
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Nu = CeReePr | 





Figure 3.2 Dimensionless representation of convection 
heat transfer measurements. 


Мо = С речева (eqn 3.5) 
Because Equation (3.5) is inferred from experimental measurements, it is termed an 


empricial correlation. However, the specific values of coefficient C and the exponents m 


and n vary with the nature of the surface geometry and the type of flow. 
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Figure 3.3 Dimensionless representation of convection 
eat transfer measurements. 


2 


IV. EXPERIMENTAL PROCEBDRE 


A. EXPERIMENTAL PROCEDURE 
To reach the thesis objectives an experimental system was established. [t is 
shown schematically in Fig. 4.1, and using water, glycol and air as a working fluids, the 
effect upon the convective heat transfer coefficient by the protuberances was 
invcstigated. 
The physical configuration of the hcat exchanger and some instrumentation 
methods will be discussed in thc following sections. 
1. Heat Exchanger 
The heat exchanger (Fig. 4.2) was formed by 14 parallel spaced plates with the 
spaces between the plates defining fluid passages. The plates are totally nested with 
each other to provide, by minimizing unsupported areas, a structurally sounder device. 
When the fluids enter the passagc, they are confronted with a unique pattern 
of interfering protuberances which come downward from the plate above and upward 
from the plate below. The turbulence created as the fluids make their sinuous path 
from inlet to outict greatly enhances the hcat transfer characteristics. Morcovcr, thc 
fact that thc rounded protuberances nest within each othcr and actually engage thc 
adjacent platc gives the heat exchanger structural rigidity which permits the use of 
thinner metallic plate material. This further enhances the heat transfer characteristics 
by lowering the metal resistancc to hcat flow. 
A brıef description of thc figures now follows : 
Fig. 4.2.a 1s an cxternal view of the heat exchanger. 
Fig. 4.2.6 1s an intcrnal view of the heat exchanger. 
Fig. 4.3 1s a detached cdgc vicw of the top plate 
configuration cmploycd in the heat exchanger. 
lig. 4.4, and Fig. 4.5 are dctached edge views of another 
type of configuration employcd in the heat exchangcr. 
Fig. 4.6 1s a detached edge view of the bottom plate 
configuration employed in the heat exchanger. 
Fig. 4.7 1s a vicw taken substantially along line a-a of 
ig. 4.2, and shows the configurations of thc platcs 


in the heat exchanger. 
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Schematic diagram of the experimental system. 


Figure 4.1 
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2. Geometric Properties Of Ihe Heat Exchanger 


PLATE 

Plate surface area : A = 

Bump surface area A Pen 

Bump cross section area : Ap, — (T r^ ®/360) — (bd/2) 
Bump base area Ay R 57/4 


= 


d 
ф/2 | е 
2 
where 
A, : plate surface area [m7] 
L : plate length | [m] 
w : plate width [m] 
Ay, : bump surface area | [m2] 
Apc: bump cross section area [m?] 
App : bump base area [m?] 
hy : bump height [m] 
b : bump base diameter [m] 
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Figure 4.2 Schematic view of the heat exchanger. 
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Figure 4.3 A detached edge view ofthe top plate configuration. 





Figure 4.4 A detached edge view of another type of configuration. 





Figure 4.5 A detached edge view of another type of configuration. 





Figure 4.6 A detached edge view of the bottom plate configuration. 
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Figure 4.7 The plate configurations in the heat exchanger. 


PASSAGE FLOW AREAS 


Hot passage flow area : a L % a 
Cold passage flow area : = у t DTA 
where 
t5 : passage height, [m] 
N, and N. are the numbers of bumps on hot and cold sides of the plate, 
respectively. 
EQUIVALENT DIAMETER 
where 
Af: passage flow area [m?] 
P : wetted perimeter [m] 
ГО PRIN USE b) 
where 
s : arc length of the bump [m] 
b : bump base diameter [m] 
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P, : hot wetted perimeter [m] 


Po : cold wetted perimeter [m] 


J. instrumentation 

Temperature measurements were made using Copper-Constantant 
thermocouples using a distilled water-ıce bath reference. A total of 16 thermocouples 
were used with individual selection being made with a thermocouple selector switch. 
This switch and an individual connector plug were installed in a tightly closed, 
insulated metal box to keep all the connecting junctions at a uniform temperature and 
thus reduce errors from connecting junctions. Four thermocouples were used on each 
inlet and outlet section of the heat exchanger. 

Pressure drops across the heat exchanger were measured using two mercury 
manometers. | 

Other instrumentation included digital flow-meters used to measure and 


control mass flow-rate so that a steady-state condition was assured. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The effects on the convective heat transfer coefficient of protuberances has been 
investigated in the fully-developed, laminar flow regime. 

The effects of protuberances are beneficial as they : 


]. Increase the structural rigidity which permits the use of thinner plate material 
and further enhances the heat transfer characteristics. 


2. Create turbulence which augments the heat transfer characteristics. 


The heat transfer correlations obtained are as follows : 
Nu = 0.113 Re® prl/3 (y/u y 
j = 0.109 Re^ 75 (yp, ) 1^ 


f = 0.549 Re "P (up) 14 


В. RECOMMENDATIONS 

The present study has been associated with the laminar flow regime. Additional 
testing should be carried out in the turbulent flow regime. | 

А data acquisition system should be employed in future testing to improve 


accuracy. 
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APPENDIX A 


SAMPLE CALCULATIONS 


CONSTANT VALUES 


Heat exchanger overall dimensions : 


length 
width 
height 
Plate length 
Plate width 
Plate area 
Bump surface area 
Bump cross sectional area 
Dump base arca 
Effective heat transfer area 
(per plate) 
Total effective heat 
transfer area (onc sidc) 
Hot passage flow area 
Cold passage [low area 
Hot side equivalent diameter 
Cold side equivalent diameter 
Flow passage height 
Bump arc length 
Bump arc radius 
Bump base diameter 
Hot wetted perimeter 
Cold wetted perimeter 
Bump numbers on the hot side 


Bump numbers on the cold side 
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209.55mm 
205.58mm 
63.5mm 
108.74mm 
104.78mm 
11393.8mm* 
108.7mm? 
226mm? 
77.1mm2 
14243.9mm2 


199414.2mm2 


141.9mm? 
106.7mm? 
2.36mum. 
520000 

3. 8559100: 
12.426mm. 
5.45mm. 
9.906mm. 
240.182mm. 
234: 79 2 пати“ 
9 

10 


(8.25іп) 
(8.094in) 
(2751m 
(4.28125in) 
(4.125in) 
(17.66in2) 
(0.16855in2) 
(0.3503in2) 
(0.1195in2) 
(27.078in?) 


(309.092in?) 


(0.21989in?) 
(0.1653in?) 
(6.093in) 
(0.07155in) 
(0.125in) 
(0.489in) 
(0.215in) 
(0.39in) 
(9.456іп) 
(9.242in) 


2. SAMPLE CALCULATIONS 
In this thesis, the following calculation procedure was uscd to calculate heat 


transfer characteristics and to obtain the correlations. 
a. CALCULATION PROCEDURE 
1. Determine the bulk temperatures, T; and t., 


ВЕКЕ) 2 
апа 


t= (t * t)/2 


2. Determine u, Ср р апі k and the Capacity rates, C, and Cy, 


Cc = mcCpc 


and 
Ch = трСоћ 


ED enne the exchanger eflectiveness, £, 


ЕЕ EL e. ul. vo) if с 
or 


DEEST AU СТ, m if en en 


4. Determine the number of transfer units, N 


Dus from EN. relations, both fluids 
Шшітікес, ( КеГб.) 


c= 1 — ехр у С (М, Чер(- СМ) 9-11 


where C, is the capacity rate ratio, C C 


min! max 


5. The overall heat transfer coefficient, U, is formed from the number of transfer units. 


UNO РА 


Cum] 


where A is the total effective heat transfer area. 
С Phe comvective heat transfer coeflicient, Aus will be 
О Ip) ra 7 


Where hg is the hot fluid heat transfer coefficient 


ће is the cold fluid heat transfer coefficient 


Di 


бу; is the wall resistance 


7. The fluid velocity, V, 1s for either fluid 
У Врт 
where Aris the passage flow area. 
8. The wall temperature, T,,, is computed from 
Ty = t+ [hp/(hy +h UT - 9 
where T is the average temperature of the hot fluid. 
tis the average temperature of the cold fluid. 
10. Form the Prandtl number, 
сым k 
the Reynolds number, 


pVD/p 


and the viscosity correction factor, 


II 


with H calculated at wall temperature. 
11. The friction factor, f, is determined from 
г = (Ау јаз 28.AT/G’v )+(1- 0% K Avy/v 7 (K ¿+1 o) 2(у И 


where v, : fluid specific volume at mean condition. 
: fluid specific volume at inlet condition. 

: fluid specific volume at outlet condition. 
G : mass velocity. 

A : total heat transfer area. 

Ar: free flow area. 

с : ratio of free-flow arca to frontal area. 
бу = 0.028 for hot side. (Ref. 6.) 

6, = 0.0178 for cold side. (Ref. 6.) 


K. : contraction loss coefficient. 


= ].3 both hot and cold side. (Ref. 6.) 


; expansion loss coefficient. 


ка 
Ке 
Кон = 0.92 for hot side. (Ref. 6.) 
Kec = 0.97 for cold side. (Ref. 6.) 


b. HOT WATER-COLD WATER STUD Y 

The experimental data ıs tabulated in Table 1. 

The experiment results are tabulated in Tables 2, 3 and 4. 
c. HOT GLYCOL-COLD WATER STUDY 

The experimental data ıs tabulated ın Table 5. 

The experimental results are tabulated in Tables 6, 7 and 8. 
d. HOT WATER-COLD AIR STUDY 

The experimental data 15 tabulated in Table 9. 

The experimental results are tabulated in Tables 10 and 11. 
e. CALCULATION OF CORRELATIONS 


ПРИТОКА ОТОН И IO OR THE NUSSELT NUMBER : 
All of the experimental data for a correlation for the Nusselt number are 
plotted in Figure A.1. Observe that the results show very little scatter, because of this, 


one may take runs 3 and $ in the water-water study and begin with a proposed form : 
Nu 7 C Ref? pr! (yg, ) !* | 

Then, it 1s only a matter of algebra to show that 
LogNu = LogC + mLogRe + nLogPr + O.ldLog(H/1t,,) 

and from the Table 4, for runs 5 and 8 respectively 


Log(10.2) = LosC + mLog(631.8) + (1/3)Log(4.79) + 0.14Log(0.892) 
Log(11.4) = LogC + mLog(756.4) + (1/3)Log(4.80) + 0.14Log(0.887) 


from which one may solve for C and m, 


C = 0.113 
m = 0.62 


> 


Hence, the correlation 1s 


Nu = 0.113 Re‘ pr1/3 ya," 


2. THE CORRELATION EOR TEIE СОРТ 
All of the experimental data for a correlation for the Colburn factor are 
plotted in Figure A.2. Observe that the results show very little scatter. Because of this, 


one may take runs 5 and 15 in the glycol-water study and begin with a proposed form : 


ј= Свет буду! апа 
Logj = LogC + mLogRe + 0.14Log(H/H,,) 


and from the Table 8, for runs 5 and 15 respectively 


Log(0.0258) = LogC + mLog(42.77) + 0.14Log(0.651) 
Log(0.0200) = LogC + mLog(83.15) + 0.14Log(0.652) 


from which one may solve for C and m, 


C = 0.109 
m = — 0.38 


ейге the sonnelation 15 


j = St Pr?/3 (n/p,,) 


ј = 0.109 Ве (p/p, y? 


where St is thc Stanton numbcr, St = h/pC,V. 


3. THE CORRELATIONEFORTTHERINIIETT OR Ро ОК" 
All of the experimental data for a correlation for the friction factor are plotted 
in Figure A.3. Observe that the results show very little scatter. Because of thıs, one 


may take runs 2 and 18 in the water-water study and begin with a proposed form : 


= С Ве? (р/р) and 
Logf = LogC + mlLogRe — 0.14Log(u/u,,) 


and from the Table 11, for runs 2 and 18 respectively 
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Log(0.051) = LogC + mLog(378.4) — 0.14Log(0.896) 
Log(0.031) = LogC + mLog(1301.0) — 0.14Log(0.897) 


from which one may solve for C and m, 


G 


0.549 
— 0.40 


Hence, the correlation is 


fz 0.549 Re" (up, ) ^-^ 
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TABLE 1 | 
THE EXPERIMENTAL DATA FOR THE WATER-WATER STUDY 


RUN AP Mh m. HOT WATER COLD WATER 
NO InH,O Т, 5 t t, 

| 1.36 0089 0.089 37.38 33.06 21.90 2022 
2 2.04 0.114 0.114 37.40 33.42 22.827 2680 
3 Su 0.146 0.146 37 50 33:501 22.94 26.53 
4 4.08 Orr Ts OS 37.28 25110 2297 2655 
5 4,76 0.190 0.190 30025 Do 2270 22 
6 5.44 0.206 0.206 о 33.66 22.55 26.04 
7 by | 0.216 / 20.216 3720 33.75 2259 — 25899 
8 6.12 02297 202223 37.10 SOT 2212 25 
9 7.48 0.248 0.248 36.90 oe 21.88 2 
10 8.16 0.2678 v 0.28 36.80 33.48 21.85 2% 
11 9.65 0.286 0.286 36.10 30:90 21.70 24.84 
12 10.47 0305 0.305 35.90 Do 2141 24.54 
13 12.24 0.330 0.330 25950) 32.44 21.23 24% 
14 13.32 0.349 0.349 35 50 32.30 21.15 24.15 
15 13.60 0.368 0.368 35.00 32.5 21.18 24.06 
16 15.64 0.381 0.381 34.60 31.80 20.877237 
17 16.32 0.400 0.400 34.50 se 20.88 2725 
18 17.68 0.412 0.412 34.35 31.68 21.00 25.67 
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АННЕ" 2 
I Ps P SPP S'OPFTIIE HOT WATER 


RUN Thh p СЪ OS 

| 35.22 1006.30 4178 721.344 
2 35.41 1006.40 4178 718.532 
сы 1006.40 4178 717.126 
4 35.50 1006.40 4178 717.200 
5 35.49 1006.40 4178 717.348 
6 35.41 1006.36 4178 718.606 
7 35.48 1006.40 4178 717.570 
8 35.39. 1006.35 4178 718.902 
9 35.21 1006.30 4178 721.492 
Ш UU 1006.26 4178 722.528 
Il 34.53 1006.01 4178 731.556 
2 3434 1005.90 4178 734.442 
Eme ur 1005.80 4178 739.844 
14 33.80 1005.70 4178 742.360 
5 8856 1005.60 4178 745.875 
№ 3220 1005.50 4178 751.240 
КОСУ 1005.40 4178 752.350 
18 33.02 1005.40 4178 753.980 


en 


k 10° 


022.152 
625.456 
625.608 
625.600 
625.584 
625.448 
023,300 
025.416 
625.136 
625.024 
624.048 
270 
025.152 
622.580 
022.500 
621.920 
621.500 
621.624 


97 
C 


— О со —) CN CA 415 (59 N 72 
— Ха 
Oe 


јр 
о -- 


SOME RESULTS OF THE WATER-WATER STUDY 


0.279 
0.273 
0250 
0.247 
0.242 
0.239 
(08755 
(02229 
09725 
0272 
0.218 
0.216 
0.214 
0.212 


0.208 


0.203 
0.202 
0.200 


T 
Ntu 


0.396 
0.364 
ше 
0.334 
0.324 
0.31) 
0.308 
0.301 
0.294 
0.289 
0.282 
0.279 
0 
0 
0.266 
0297 
0756 
0:290 


ТАВИ Е 3 


U 


737.40 

919.48 

1057525 
1199.68 
1293.06 
137945 
1593 03 
1441.46 
1525.18 
J 615. Pe 
1687.52 
1780.73 
1903.12 
1988.80 
2050-55 
2049.47 
2143.50 
2185.54 
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hh 

1510.2 
1894.2 
2145.0 
2494.2 
2696.7 
2884.4 
2914.9 
3021.0 
3205.3 
3400.2 
3565.8 
3774.5 
4050.6 
4245.3 
4386.9 
4425.1 
4599.6 
4696.6 


hc 


15102 
1894.2 
2145.0 
2494.2 
2096.7 
2884.4 
2914.9 
3021.0 
220575 
3400.2 
OO 
IES 
4050.6 
4245.3 
4380.9 
4425.1 
45926 
4696.6 


29.64 
30.11 
30.12 
30.08 
2993 
29 
29.60 
29501 
29009 
29.32 
28.90 
28290 
28.30 


28.23 


28.10 
27 
2707 
27.68 


LA 
@ 
a 


u EI ER 


на | 
Soa 


15 


Pr 


4.82 
4.80 
4.79 
4.79 
4.79 
4.80 
4.79 
4.50 
4.52 
4,90 
4,92 
4,96 
4.95 
5.01 
08 
5.06 
24017 
EID 


St 


0.0040 
0.0039 
0.0035 
0.0035 
0.0034 
0.0033 
0.0032 
0.0031 
0.0031 
0.0050 
0.0050 
0.0029 
0.0029 
0.0029 
0.0028 
0.0028 
0.0027 
0.0027 
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| TABLE 4 
SOME RESULTS OF THE WATER-WATER STUDY 


Re 


293.3 
378.4 
484.5 
568.7 
631.8 
683.2 
715.7 
756.4 
816.5 
378.0 
928.9 
986.8 
1061:1 
1118.3 
1173.8 
1205.4 
1262.9 
1301.0 


i 


0.891 
0.896 
0.896 
0.894 
0.892 
0.892 
0.889 
0.887 
0.886 
0.886 
0.890 
0.889 
0.890 
0.891 
0.892 
0.891 
0.894 
0.897 


0:055 
0.051 
0.049 
0.048 
0.043 
0.042 
0.041 
0.040 
0.040 
AVES 
0.040 
0.038 
0.058 
0.037 
0.036 
0.034 
0.033 
0051 


| 

0.0115 
0.0112 
0.0099 
0.0098 
0.0095 
0.0094 
0.0091 
0.0089 
0.0088 
0.0087 
0.0086 
0.0085 
0.0085 
0.0084 
0.0083 
0.0081 
0.0081 
0.0080 
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| TABLE 5 
THE EXPERIMENTAL DATA FOR THE GLYCOL-WATER STUDY 


RUNAP my m. HOT GLYCOL COLD WATER 
МО inHg T, il t U: 

| 0.3 0.098 0.089 :36.5 ие т 210 22.8 
2 04 0.125 0.114 36.5 33.80 21000 
3 0.5 0160 0.146 36.5 33.97 21.0 255 
4 06 018% 017 364 33.93 21.0 226 
5 08 0.209 0.190 364 33.99 21.0 Ж 
6 09 0227 0206 364 34.03 20, 22.4 
7 10 0237 0216 364 34.09 20.9 22.4 
8 11 0.251 0.229 365 34.24 20.9 22.4 
9 13 0272 0.248 36.5 34.30 21.0 22.4 
0 14 09.293 0.267 36.5 34.33 210 224 
и 16 0.314 0.286 365 34.37 210 224 
ІМ 19 Оо 36.5 34.44 21.0 Ж 
и ооо 36.5 3447 21.0 “бәй 
4 23 0383 0.349 36.5 34.50 210 223 
5 26 0.404 0.368 36.5 BASS О ame 
6 27 0.418 0.381 364 34.46 209 221 
7 29 0.439 0.400 36.4 34.50 209 22.1 
18 31 0453 0412 364 34.53 209 224 
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LA 
QE 
Z 


о CO мм Сс tA A (59 N — 


Pm И и = 
СО -3 QDR л HB O кд = © 


= 
e 
SI 


БӘРІ? 
55.15 
35.24 
BS 2 
25.20 
3927 
59-29 
БЕЗІ 
35.40 
35.42 


35.44 


35.47 
35.49 
33:810) 
Br 
35.43 
35.45 
35.47 


p 


1109.94 
1104.9] 
1104.96 
1104.94 
1104.92 
1105.05 
1 105992 
1105.01 
110509 
1105.05 
1105.02 
1104.96 
1105.03 
1105:05 
1105.05 
1104.96 
1105.00 
1104.96 


TABEE 6 
БОБОЙТ НОТ ОД УСОТ, 


2452 
2452 
2452 
2452 
2452 
2452 
2452 
2452 
2452 
2452 
2452 
2452 
2452 
2452 
2452 
2452 
2452 
2452 
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it 10% 


1.164 
СІБЕ 
ШЕ 
1.162 
наше 
MISS 
165 
152 
12156 
i. 
| 
1 
| 
1 
l 
1 
l 
1 


149 


.148 
.147 
.146 
. 145 
.144 
.149 
.148 
.147 


k 10° 


254.5 
254.5 
254.5 
254.5 
254.5 
254.5 
254.5 
254.5 
254.5 
254.5 
254.5 
254.5 
254.5 
254.5 
254.5 
254.5 
254.5 
254.5 


Z 
N 


Oo CO — C C d. о YR 
СУС; 


к. ке шч =— =— 
RW < 


ЗОМЕ. ВЕЗЫУГТ ОЕ ТИВ ССОО» 


ШІ? 
0.174 
0.163 
0.160 
0.156 
0.153 
0.149 
0.145 
0.142 
0.140 
0957 
0153 
ШІ | 
Zo 
ШІСІ 
OTIS 
25 
МЕЛІ 


"tu 


022219 
0.204 
0.189 
0.185 
0.180 
0.176 
0.170 
0.165 
0.161 


0.159 


OSs 
0.150 
0.147 
0.145 
0.142 
0.140 
ШӘП 
Оо 


ТЕ 


U 


270.01 
314.73 
372.64 
426810 
462.88 
490.38 
495.55 
90925 
538.30 
EAS 
39 7:9 
617.24 
655.32 
683.71 
706.09 
720.08 
TI 
15228 
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hy 


332.0 
382.02 
457.53 
525.49 
568.88 
602.11 
608.56 
624.59 
660.49 
7098 


To dais 


153958 
801.67 
836.57 


864.60 


884.13 
907.08 
92158 


ой 
1894.2 
2145.0 
2494.2 
DOUO T 
2584.4 
2914.9 
202177 
pos) 
3400.2 
Бобо о 
3774.5 
4050.6 
4245.3 
4586.9 
4425.1 
4599.6 
4696.6 


W 

24.77 
24.74 
25.00 
24.95 
24.86 
24.83 
24.98 
24.99 
24.99 
24.89 
24.86 
24.81 
24.66 
24.65 
24.72 
24.74 
24.60 
24.55 


Z 
Z 
E 


= Ono у су сл с ко == 
c “ж 
Qe 


pi — 4 
м = 
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Pr 


112.18 
[1201 
КУ 
1225 
112.11 
111.54 
111.37 
111.68 
MESS 
152 
111.42 
111.68 
111.61 
111.54 


111.49 


Ши 59 
52 
DIES 


TABLE Š 


St 


0.00130 
0.00125 
0.00116 
0.00114 
0.00111 
0.00108 
0.00104 
0.00101 
0.00098 
0.00098 
0.00095 
0.00092 
0.00089 
0.00088 
0.00087 
0.00086 


0.00084 


0.00083 
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Re 


97 
PICS, 
DD 
36.45 
42.77 
46.57 
48.80 
OZ 
35.86 
60.19 
64.55 
68.70 
74.47 
78.82 
55715 
Soy 
90.05 
I2 2 


ШИ, 


0.650 
0.651 
ШО 59 
0.654 
0.651 
0.001 
(A 
08699 
DEOS 
0.656 
0.656 
0:055 
0.649 


0.650 


0.052 
0.656 
0.647 
0.647 


0.142 
0.115 
0.090 
0.090 
0.094 
0.090 
0.090 
0.090 
0.087 
0.084 
0.082 
0.081 
0.080 
0.080 
0.080 
0.079 
0.077 
0.077 


0.0302 
0.0290 
0.0270 
0.0266 
0.0258 
0.0249 
0.0241 
0.0234 
0.0228 
0.0226 
0.0220 
0.0213 
0.0207 
0.0205 
0.0200 
0.0200 
0.0195 
0.0193 


ТАВ 9 
THE EXPERIMENTAL DATA FOR THE WATER-AIR STUDY 


mp HOT WATER COLD AIR 
| 2 ti L, 
36.8 36.8 19.6 36.5 


36.8 80.7 19.6 36.3 
36.3 86.7 1926 36.0 
30.) 36.8 19.6 36.9 
pon 36.8 19 DN 
36.9 36.5 19.8 3846 
S 36.7 19,8 2959 
36.9 ЕСУ ша М2. 
O70 36.8 2 3589 
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0.983 
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0.953 
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0:956 
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0.918 
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0.901 


E 
E 


to, сз > 
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00 


2.78 
20l 
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U 


39.74 
31:55 
з о 
69.66 
80.72 
88.31 
97.86 
104.47 
113.69 
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Dh 

2696.7 
2696.7 
2696.7 
2696.7 
2696.7 
2696.7 
2696.7 
2696.7 
2696.7 


40.58 
52.99 
61.23 
"ION 
83.43 
91.56 
101.88 
109.06 
119.14 
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Figure A.3 A correlation for the friction factor. 
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APPENDIX B 
DIMENSIONAL ANALYSIS 


The dependent variable for the process is the convection heat transfer coefficient, 
h. For incompressible low-speed flow the independent variables that determine the 
heat transfer coefficient are the fluid velocity, V, a linear dimension (e.g., the equivalent 
diameter, D ), and the fluid properties of thermal conductivity, k, viscosity, H, specific 
heat, Ср and density p 

The independent dimensional quantities to be used in the analysis are mass, М; 
length, L; time, ©; and temperature, T. The variables, their symbols, and their 


dimensional equations are listed ın Table 12. 


TABELLE 
VARIABLES FOR CONVECTION 


VARIABLE SYMBOL DIMENSIONAL EQUATION 
Equivalent diameter D [L] 
Thermal T k [ML/O?T] 


of the flu 

Velocity of the fluid  : [L/O] 
Density of the fluid [M/L?] 
Viscosity of the fluid [M/LO] 


Specific heat at [L?/O*T] 
constant pressure 


Heat transfer [M/O?T] 
cocfficient 





There are seven physical quantities and four primary dimensions. We therefore expect 
that three dimensionless groups will be required to correlate experimental data. To 
determine these dimensionless groups, Write as a product of each of the variables 


raised to an unknown power: 
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П = га ко ус pd ye ES hS 
and then substitute the dimensional formulas from the tabulation above: 
П- памуео тема? 9  млер2? от мезпе 


For Il to be dimensionless, the exponents of each primary dimension must separately 
add to zero. Equating the sum of the exponents of each prunary dimension to zero 


yields the following set of equations: 


ше Ре go for M 
a + b + c — 3d — e + 2f = 0 for L 
ӘН СЕ БУЕ 36 = 0 for Q 
ES fumo) for T 


Since there are seven unknowns but only four equations, three of the ролен ере 
arbitrarily selected for each dimensionless group. The only restriction 1s that each of 
the arbitrarily selected exponents be independent of the other. This requirement IS met 
if the determinant formed by the coefficients of the remaining terms does not equal 
Zero. 

Since the convection heat transfer coefficient is the dependent variable that we 
eventually want to evaluate, set its exponent g equal to unity in the first dimensionless 
group, Il, To simplify the algebraic manipulations, arbitrarily set c = d = 0. This 
means that in the first dimensionless group velocity and density will not appear. 

Solving the equations above simultaneously, we obtain а = 1], b » — l, ande 


= f = Q. The first dimensionless group is therefore 
П, = hD/k 


which we recognize as the Nusselt number. 
For the second dimensionless group, we set g = 0 to avoid having the 
dependent variable h appear again. Arbitrarily, we leta = I and f = 0. Simultaneous 


solution of the equations above with these choices yields the second parameter, 
Il, = VDp/u 


which we recognize as the Reynolds number, Re. 
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For the third dimensionless group, we wish to have the specific heat, which has 
not appeared in any previous group, to be included and therefore set f = 1. Similarly, 
we do not wish to have the heat transfer coefficient or the diameter appear again and 


therefore let a and g equal to 0. This yields the third dimensionless group, 
Il, = Сым 


which we recognize as the Prandtl number, Pr. Thus the heat transfer coefficient in T, 


can be functionally related to the dimensionless groups II, and Il, in the form 


Nu = f Re, Pr) 
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APPENDIX C 
UNCERTAINTY ANALYSIS 


The uncertainties for the major variable in the experiments were calculated in 


accordance with the method described by S. Kline and F. McClintoch (Ref. 9.). The 


estimates of the uncertainties were made conservatively. As a result, there is 


considerable confidence in the uncertainties as calculated. 


The following equations were used to calculate the uncertainties. 


ЈЕ 


G3 


Š. 


9: 


дтутъ - ГОт/т))? + (дтутЬ 12 


, деј Су = [(0m/m)? + (дС С) 1? 


де = L (OCy/ CS? -- (OTT, + (ат /т,) + (дол, 12 


. NA NA, [oG G + (деје)2 112 


. 8UJU = [KON у Хур 2 + (CH/CH)? + (A/A) 1? 


V/V = [ (öm/m)? + (öp/p .)? + (ВАНА? 1? 
дн = [ (20/0)? + (Gr, Jr. y? Y 
üPr/Pr = (С/С)? + (дни)? + (дкло? 192 


ORe/Re = [(0p/p? + (0V/V? + (0D/D? + (Opu) JP 


10. 9Nu/Nu = [ (móRe/Re? + (noPr/Pry? ]1/2 


11. Əf/f 


12. ijj 


| 


(möRe/Re) 


| 


(тд Ке/ Ке) 


The following quantities had their uncertainty calculated by dividing the estimated 


error in the quantity by the value of the quantity and these uncertainties are considered 


constant for the range of values in this experiment. 


ШЕ 


Quantity Uncertainty 


0.003 
(005 
0.006 
0.005 
0.003 
0.002 
0.008 
E 0.005 


Der 0.003 


зсоғсЕ О > > 
| =; 


+ 


The following quantities had the worst case error 


Quantity Error 
Т, И 
T, + O 
t +0? 
t. zu 
AT әс 


Then, the uncertainties for the following quantities were calculated using the equations 


above, 

Quantity Uncertainty 
Ty _ 0.008 
Ch 0.01 
с 0.012 
Ns 0.015 
U 0.018 
h 0.018 
V 0.01 
Pr 0.009 
Re 0.013 
Nu 0.009 
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